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Abstract. The aim of the paper is to relate computational and arithmetic questions 
about Euler's constant 7 with properties of the values of the q-logarithm function, with 
natural choice of q. By these means, we generalize a classical formula for 7 due to Ra- 
manujan, together with Vacca's and Gosper's series for 7, as well as deduce irrationality 
criteria and tests and new asymptotic formulas for computing Euler's constant. The 
main tools are Euler-type integrals and hypergeometric series. 



We recall the definition of Euler's constant: 



lim f logn ) 



7 

^ 0.57721566490153286060651209008240243104215933593992 . . . 

and the fact that the (expected) irrationality of 7 has not yet been proved. 

Recently, the first author gave [9], [10] a construction of Z-linear forms involving 7 
and logarithms; namely, he proved that 

rfan^n G Z + Z7 + Zlog(n + 1) + Zlog(?i + 2) H + Zlog(2n), 

where 

//^ll^^rVil^dxdy 

{\ - xy)\\o<gxy\ 



[0,1] 



2 



^ /" ('t(t + l)---(t + n)') 
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and dn denotes the least common multiple of the numbers 1, 2, . . . ,n. This type of 
approximation allows deducing several (ir)rationality criteria for Euler's constant [9] . 
In a sense, what follows is inspired by considerations in [9] . 
A special function that we will require is the g-logarithm 

(1) ln,(l + .) = ^ ^ , \z\<q, 

1/=! 

for q real with \q\ > 1 (in fact, our main choices will be = 2 in Sections 2-4 and 
q = 3 in Section 8, but we let > 1 be any integer in Section 9). The function (1) is 
a qf-extension of the ordinary logarithm function 

since 

lim((Q-l)lnq(l + 2;)) =log(l + ;2), \z\ < 1. 

q—>l 

We also mention that the irrationality of the values of the g-logarithm (1) for q G Z\{0, 
±1} and 2; e Q, 2; ^ 1, is known [3]. 
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1. Gosper's acceleration of Vacca's series 
The following expression for 7 is known as Vacca's series [11]: 

(2) , = 

n=l 

where log^a; = (loga;)/(log(j) is the ordinary logarithm base q. Gosper [6] has trans- 
formed this series into a more rapidly converging series of positive rationals. His 
method (see Section 6) is to use partial summation to write (2) as the double series 
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apply Euler's transformation to each inner sum on k, obtaining the double series (7) 
below, then triangularize it. In this way, he arrives at a version of series (4). Here 
we give a different proof of (4). To begin, we reverse the order of summation in (3) 
and write the resulting double series as the integral (5), which is a form of an integral 
that Bromwich attributes to Catalan ([4], p. 526). 

Theorem 1 (Gosper's Series). We have the series for Euler's constant 



i/=2 k=l ^ 

Proof. After expanding 1/(1 + x) in a geometric series, termwise integration shows 
that the formula 



(5) 



is equivalent to (3) with the order of summation reversed, which is easily justified by 
grouping terms in pairs. Integrating by parts K + 1 times gives 

^'h^'^^^'^Tk')^ 'o (l + a:)^+2Z.2^(2. + i)...(2^+K) 
Since the integral is less than 

J„ (l+x)'f+2;^2" K + l\ 2'f+V' 

as K tends to infinity (6) becomes 

oo oo ^ 

^ nv+k+1 a^+k^ ■ 

k=0 v=l ^ \ k ) 

The inner series with = sums to 1/2, and in the series with A; > we may collect 
terms (triangularize) as in (4), completing the proof. □ 

Remark. We can regard the series on the right of (2), (4), (7) as base q = 2 series 
for Euler's constant. Gosper indicated [6] the following extensions of (2), (7) to base 
q = 3: 



oo I 1 I 2nn oo 7r(fc — 1) 

(8) ^^2^ Llog3nhcos— ^^^^ cos^ 



6 



n A^A^^u+{k+i)/2ri''+k\ 

v^l k=0 ^ \ k J 



(see Sections 5 and 6 for proofs and generalizations), which give faster convergence 
than formulas (2), (7), respectively. 
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2. A g-LOGARITHM APPROACH 

Lemma 1. Let q be real, q > 1. Then for all integers n ^ and k, with < k < q^~^^ , 
the following identity holds: 



V ^ =-ln (l -) 



(9) 

Remark. For k — 0, the series in (9) equals 

1 q-in+i) 



- 1 _ q-(^+^) _ 1 

Proof. We have 



^ 1 ^ 1 1 ^ 1 (-/c)^-i 

^ q"" + k~ ^ g^l + fc - o-'" ~ 2-^ qy 2-^ q^ki^-^) 



oo ^ oo 



/i=l i/=n+l /i=l 

A* — 1 



as required. □ 



Remark. P. Sebah [8] has pointed out that in order to compute the qf-logarithm, the 
following formula is much more efficient than (1) or (9): 

\n,{l + z) = zJ2-r^+rN{z), where rjv(2:) = ^ 

with N being any positive integer. The proof of the formula is similar to that of 
Lemma 1. 

Let m be a non- negative integer. Using the rational function 



for each ^ define the convergent series 

oo 

(11) 5.,„ = 5^(-l)*i?„(2^ + 0. 

t=0 
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Lemma 2. We have 



(12) 



k=0 .1 = 1 ^ ^ fe=0 



(We use the standard convention that any sum Ylj=a ^^^^ if a > 6.) 
Proof. Indeed, 

_ - - (-1)^+^(7) _^fm\^ 



i=0 ^"'^ A;=j j=0 / ^fc=0 fe=0^ 

and (12) follows. □ 

Lemma 3. For the series Sum, we have the representations 



(-1)' 



(13) 



1 + 



Proof. Replace t by t — 2^ in (11) to get the first equality in (13). For the second, 
use the Binomial Theorem to write 

Rm{t) = [ {l-x)'^X*-^dx. 

Jo 

Substitute this in the series for Si^^m, interchange summation and integration, and 
sum the series 



(14) 



E(-i) 



X 



2''-l 



1 + X 



arriving at the second equality in (13). To justify termwise integration, replace the 
geometric series (14) with a finite sum plus remainder rN{x) = {—l)^x^'^'^"~^/ 
(1 + x), and note that 



rjv(x) dx 



< 



1 



as iV ^ oo. □ 



Our final object is the double series 



(15) 



'■n,7n — / ^ >~'i>,m-i 
u=n+l 
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where n ^ 0. By Lemma 2 and formula (3), for the quantity (15) we have 

/ °° " \ (-1)'' °° (-1)'' 

v=V k=0 j=l ^"^ ^ fe=Oi^=n+l 

=2'»,-2-5:(i:-i:)ui-e(;) e ^ 

1^=1^2=1 /=1 ^ j = l ^''^ i/=n+l 

For simphcity, assume now that m ^ 2^+^ Then using Lemma 1, 



n 2"-! 



(16) /,,^ = 2™-(7 + nlog2)-2-5^5^ ^ ' 



I 

u=l 1=1 



We can summarize the result of these computations in the following way. 
Lemma 4. Ifm^ 2*^+^, then 



where 



and 



satisfy 



K,^ = 2-nlog2 + g (^) g ^ m. (l + I, 



T 



i/=l 2=1 



d„L„^ eZlog2 + Zln2( 1 + ;^ ) +Zln2ri + |-) + • • • + Zln2 ( 1 + ^ 



2n J \ 2" y \ 2" 

and (i2"^n,m £ ^• 

Lemma 5. For the double series In,m) have the integral 

(18) In,m= /'iiZ^ g x^-Mo;, 

which generalizes Catalan's integral (5) forj. 

Proof. Integrating termwise (since the series in the integrand is dominated by a 
geometric series, the justification in the proof of Lemma 3 applies here) and us- 
ing Lemma 3, the formula follows from the definition of /n,m- The case /o,o = 7 
(see (16)) is indeed the integral (5). □ 
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Lemma 6. For m a positive integer and r ^ 1 a real number, we have the bounds 
(19) 

and, for < a; < 1, 

(20) x'^'^-^l - x)"^ < 4 



(r + 1) 



r+l 



Proof. Euler's beta integral 



[\'^-\l-xf-^dx^^^^^^^^, Rea>0, Re/?>0, 
Jo r(a + /3) ' 

gives 

x^"^-Hi - xr dx = =±Y\(i + r^y\ 

Jo ^ ' r(rm + m+l) rm j J ' 

denote the product by n(m, r). Using r ^ 1 to get 



J + rm m \ r 

we obtain 



n(m + l,r) 1 -TT i + rm If 1^ ^ 

- ^n-^^ —>—rr(l + - 



n(m,r) r + 1 j-J-^ j + r{m + 1) r + l\ rj (r + 1)^+-'- ' 
and by induction on m it follows that 



n(m, r) > 



(r + 1)^+1 



since the inequaUty holds for m = 1. This proves (19). 
If < a; < 1, then 



(l-x)'" < (x'^(l-x))"' ^ 
(r + 1)1+1/^ 



(r + 1)'"+! 



m—l/r 



(r+l) 



r+l 



and (20) follows, using r ^ 1. □ 
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Lemma 7. If n ^ is an integer and r is a real number satisfying 1 ^ r ^ 2^^^+^, 
then we have the bounds 
(21) 



^ \ m+l / J. 



2r(m+l) V(r + 1)'"+V "-^ - " \^(^ + i)r+i 



T \ IT 

< In,m < 6 • , , 1 for 



2n+l 



In particular, for r ^ 1 fixed and m = [2'^'^^ /r\ ^ 1, t/ie quantity In,m is positive and 
decreases exponentially at the rate 

(22) lim ^"S-^"-L^"+'/^J = rlogr - (r + 1) log(r + 1). 

Proo/. Put m = [2"+^^ Lemma 5. Then m + 1 > T'^^ jr, hence 



C (L_^a;2"+^-i dx>- f' - x)-+^ dx, 



and the lower bound in (21) follows from (19) (with m + 1 in place of m). 
Using m ^ 2"+-'^/r and relations (20) and (5), we have 



Tl^ 1 



Jo 1 + ^1 

< .—^(1 - (i + f: -^^-^) 

and the upper bound in (21) follows from the inequality 4(log2 + 7) < 6. Since (21) 
implies (22), we are done. □ 

Remark 1. If 2"+^/r is an integer, then we may replace m + 1 by m in the lower 
bound in (21). 

Remark 2. As r — 00, the right-hand side of (22) equals — logr + 0(1), hence its 
absolute value tends to infinity with r. 

3. Irrationality tests for Euler's constant 

Theorem 2 (Rationality Criterion for 7). The fractional part of (i2"-^n,m equals 
d2^In,m for all n and m with 

(23) n sufficiently large and 2^~^ ^ m ^ 2"+\ 



if and only if Euler's constant is a rational number. 
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Proof. With 2"-^ ^ m ^ 2''+^ set r = 2"+Vm, so that 1 ^ r ^ 4. By Lemma 7, 
we have /^^^ < 6p(r)"^ = 6p(r)^"'^/'", where p(r) = r'^/(r + 1)'^+^. Straightforward 
verification shows that p(r)^/'' is increasing for r > 0, and satisfies 

0<p(r)2/^<- for < r < ro = 5.6213305349 ... , 

e 

where p(ro)^/''° = 1/e. Therefore In,m < 6(p(4)i/2) and p(4)i/2 < 1/e. Since the 
Prime Number Theorem imphes that for any £ > we have (In = 0((e^"''^)^) as 
N — > oo, we conclude that < d2nln,m < 1 for n and m satisfying (23). 
Now multiply (17) by and write the result as 

(24) d2nLn^rn = d2nln,m + d2r<-{An,m — 2*^7). 

Since (i2"^n,m € Z by Lemma 4, the theorem follows. □ 

Remark. The argument shows that the Criterion holds with (23) replaced by the more 
general conditions n ^ n(ri) and [2'^''"-^/ri] ^ m ^ 2"'"''^, where ri is any number 
between 1 and tq. 

Theorem 3 (Irrationality Test for 7). // the fractional part of d2"Ln,2^ satisfies 

(25) K"^n,2"}>6-(^) 

infinitely often, then Euler's constant is irrational. In fact, the inequality for a 
given n > implies that if ^ is a rational number, then its denominator does not 
divide 2^ d2n . 

Proof. It follows from [7], Theorem 13, that d2N < 8^ for all integers N ^1 (see [9], 
Lemma 3). Combined with the upper bound in Lemma 7 for r = 2, this gives 

ry cyn — 1 rj 71 — 1 

4 \ / 128 ^ 

< d2"/n 2" < 8^ ' ' • 6 ■ ( — =6 



2"^ 



27y V729 

for n ^ 1. The first part of the theorem follows by the Rationality Criterion for 7, 
and the second part by formula (24) with m = 2". □ 

Remark 1. Condition (25) in Theorem 3 can be replaced by the simpler condition 

|w2"-^n,2" I > 5 

if n > 1. For n > 4, this follows from the inequalities 128/729 = 0.17558299 . . . < 1/5 
and 6 • (128/729)^^ < 5"^^. For n = 2, 3, 4, it follows from the proof of Theorem 3 by 
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replacing the bound < 8"^ with the exact values d4 — 12, dg = 840, diQ — 720720, 
for which ■ 6 ■ (4/27)2" < 5~^^~\ 

Remark 2. Using the Rationality Criterion for 7, we may generalize the Irrationality 
Test as follows: 7 is irrational if there exist infinitely many integers n and m such 
that 2^-^ ^ m ^ 2^+^ and 



(26) {d2nLn,m\>8 ■ 6 • { ^^_^-^y+i J > where r = 



2 



n+l 



m 



For example, the case m = 2"^ is Theorem 3, and the case m = 2^~^^ can be stated: a 
sufficient condition for irrationality of 7 is that 

|(i2"i>n,2"+i ) > 6 • 2 ■ infinitely often. 

For any e > 0, we may further refine the Test by replacing 8^" ^ • 6 with e^^'^^^'^'^ 
in (26). 

Example. P. Sebah [8] has calculated that 

{^212^12,212} = 0.178346164... . 

Since the last number exceeds 5"^ , if ^ is a rational number a/b, then b is not a 
divisor 0/ 2^°^^ • ^4096; in particular, \b\ ^ 4099 (since the numbers 4097 = 17 • 241 
and 4098 = 2 • 3 • 683 divide c/4096)- 

Remark. Similar irrationality tests for 7, but using ordinary logarithms instead of 
2-logarithms, are proved in [9], where an example shows that no divisor of 
(10000) ' '^20000 can be a denominator of 7. The present example yields additional 
information because of the high power of 2. 

4. Computing Euler's constant 
Theorem 4. The following asymptotic formulas as n — > 00 are valid: 

(27) 7= ^"^"";2:>""^ +Q(2-^-^"), 

^n,2" ~ -£'n,2" „ / / 2 



(28) 7= "V +^ 27 



^n,2"-i ~ -^n,2"-i ^( f 128 \ 



(29) - +0^^— j 

(30) 7=i^n,i-nlog2 + 0(2-). 



EULER'S CONSTANT, g-LOGARITHMS AND RAMANUJAN 



11 



Proof. For (27), put m = 2^^^+^ in (17) and use the upper bound in (21) with r — 1. 
For (28), (29), take m = 2'^,2'^~^ and r = 2,4, respectively. Finally, let m = 1 in 
Lemma 4 to get In,i — 2'y + 2nlog2 — A^^i- From (21) with r = 2"'+-'^, we obtain 

< / , < 6 ^- I - - ( 1 + 1 < 0-"+^ 

U ^ Jn,l ^ D ^2n+l _^ 2)2"+' + l ^ 2^+^ + 1 V 2"+^ / 

for n ^ 0, proving more than (30). □ 

Remark. Using (28) and his calculation of Ln,2" for n = 12, P. Sebah [8] has com- 
puted 7 correct to 4631 decimal places, which is the accuracy (2/27)^^^ predicted. 
Formula (27) gives almost 60% more digits of 7, but requires computing about twice 
as many 2-logarithms. Similarly, (29) yields nearly 39% fewer digits than (28), but 
involves only half as many 2-logarithms. Letting m continue down to 1, we reach 
formula (30), which doesn't involve any 2-logarithms. 

5. RAMANUJAN'S base q INTEGRAL FOR EULER'S CONSTANT 

For integer q ^ 2, Ramanujan (see [2] for proofs and references) gave the formula 
for Euler's constant 



(31) 



which reduces to Catalan's integral (5) when q = 2. Equation (31) implies the 
following base q generalization of Vacca's series (2) (see [2], Theorem 2.6): 



00 



(32) ^^^.nLlog^nJ 



n=l 



where 

q — I if q \ n, 



(33) an = an,q = 



-1 a q\ n. 



The case q = 3 gives the first equality in (8). 

Let e = eg be a fixed primitive q-th. root of unity (for instance, e = e^"^*/^). Using 

q-l q-l 

1=0 1=1 



we can represent formula (32) as 



q-l 00 1, ^ ^\^nl 



(34) , = ^^Llog,nJ.' 



1=1 n=l 



which can be regarded as another base q generalization of Vacca's series (2). 



12 



J. SONDOW AND W. ZUDILIN 



6. A BASE q GENERALIZATION OF GOSPER'S SeRIES 

Following Gosper's proof [6] of (7), we write (34) as 

Q — 1 oc oc 



7 



EEE 



1=1 ^=1 fc=0 + ^ 1=1 ^=1 ^ 



and apply Euler's transform (see, e.g., [1], Section 2.4, formula (1)) 

1 



7 



z = 



(1 - zY 



■ 2F1 



a,7 - /3 
7 



—z 



l-zP 



yielding 



q — l 00 



1 1 



7 



g— 1 CXD ^ CX3 /7y,7\— 1 

1=1 u=l ^ k=0 



— e 



q'' + k\ ' (-l)''r''^ 



(1 - 60 



Substituting x = 1 into the expansion 

1 -X'? 



g-l 

n(i-6%)= 



1 -a; 



= l + a; + a;^H \- x 



q-l 



leads to 
(35) 



q-l 

11(1-^') 

J = l 



(where the empty product in the case q = 2 equals 1). Therefore, we arrive at the 
following result. 

Theorem 5. We have the base q ^ 2 accelerated series 

(-i)'x.(^) 



(36) 
where 



i/=l k=0 



q-l q-l 
1 = 1 j = l 

a symmetric polynomial in the roots of the polynomial 1 + x + x'^ + • • • + x*^ , is an 
integer-valued function independent of the choice of primitive q-th root of unity e. 
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Remark. From (35) we have 

f-i 

(1 - 60 



and it follows that 

k+l / \ ^=+1 

Thus we have the relation 

9 



limsup|xg(/c)r/'' ^ 



fe^cx) 2 sin I 

which shows how the convergence of series (34) has been accelerated in (36). 

Examples. Series (36) is a generalization of Gosper's series (7), which is the case 
q = 2. When g = 3, taking e = e^^*/'^ we obtain 

_ g27rA;i/3 _ 2(A;+l)/2g7r(A;+l)i/6 _|_ g-27rfei/3 . 3(fe+l)/2g— 7r(A;+l)i/6 

= 3('=+^)/2.2cos^^^^^, 

which proves the second equality in (8). When g = 4, we take e = e'^*/^ = i and 
obtain 

XA{k) = i''2''+\l + + (-1)'=(1 - + + {-i)''2''+\l - 

4 

which is an integer- valued function of k. 

7. Base g integrals 

There are several ways to generalize integrals (18) for In,m and (31) for 7 simulta- 
neously. One way is to define 



(37) 



Then /n,m,2 = -fn,m and /o,o,g = 7) and the analog of Lemma 7 with r = q and 



m = g", 



(38) hm l^S%:^i = 5logg-(g + l)log(g + l), 
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can be easily derived by the methods of Section 2 (cf. Lemma 9 below). We may 
expand integral (37) as in the proof of Lemma 4. Namely, with Rm{t) defined in (10), 
and at in (33), we have 

q 1 (q — l)—x — x^ — --- — x^~^ + 

— = -^^ = > (ytx\ 

l-xi 1-x l-xi ^ 

t=0 

hence 

oo oo „i oo oo 

/n,m,, = J]^t E / x'^''+'-\l-xrdx= J2 + 
t=0 v=n+l •^'^ v=n-\-l t=Q 

m / \ oo oo 

j=Q / i/=n+l k=i ^ 

Unfortunately, the last representation shows that In,m,q involves not only logarithms 
and Qf-logarithms but also "generalized Euler constants" 

oo oo 

(39) 7... = EE^' i = 0,l,...,Q-l; 

u=l k=j 

the equality 70, q = 7 follows from (32) just as (3) follows from Vacca's series (2) (the 
case q = 2). We may avoid the appearance of generalized Euler constants 7^^^ 7^ 7 in 
the particular case q = S (details are in Section 8), but not for > 3 (see Section 10). 

However, there is a way to generalize integral (18) for In,m that solves this problem, 
namely, by defining 

(40) .;„„ = /"(^-^)(,-i^)'" f .«-d. 

= r (g-i-^-^' ^'-"r^' V x^'-'dx 

with the properties /4,m,2 — ^ri,m and /o,o,g = 7' (37). For integral In^rn,q in (40), 

we may generalize the results of Sections 2 and 4, but the poor asymptotics of the 
integral make these generalizations useless for irrationality testing. For example, the 
choice m — {q^ — l)/{q — l) gives the inclusions 

dcir^In,m,a e Z + Z7 + Z log 9 + Z In, (^1 + + Z In, (^1 + + • • • 

/ - 2 

+ ZlnJl + ^ 



(see Lemma 10 below) and the asymptotics 



(see Lemma 13); this makes it impossible to use In,m,q Si'^® irrationality tests for 
g > 2 similar to those of Theorems 2 and 3 for g = 2. Instead we use I'n,m,q extend 
Theorem 4 to a base g > 2 asymptotic formula for 7; details are in Section 9. 
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8. Base 3 irrationality testing 
Taking g = 3 in (31) gives one of Ramanujan's formulas for Euler's constant 

(see [2], Corollary 2.3 and the equation following (2.7)), and taking g = 3 in (37) gives 
the base 3 integral 

(42) ,_3 = ^'_J±^(1_,,™ f 

" i/=n+l 

It turns out that, for m a multiple of 6, we can extend the results of Sections 2-4 to 
base g = 3, as follows. 
Since 

}'^-'^^\ =28-34x + 21x^-7x^ + x^- ^ 27 

1 + X + X'' 1 + X + x'' 

by induction we have 

(43) ,^i^=Q.(.). 



1 + X + x"^ 1 + X + x'^^ 

where Qk{x) is the polynomial defined by the recursion 

Qo{x) = 0, 

(44) 

^ ' Qk+i{x) = {l-xf Qk{x) + {-2,f^ {28- Mx + 21x'^ -Ix'^ + x^) for A; ^ 0. 
From (42), (43), (41), for m = 6A; we obtain 

\x^'^-^dx 



(45) /„,^,3 = j\2 + X) {q,{x) + yH^) (E - t.} 



„1 oo 

+ / {2 + x)Qk{x) J2 



x^"-^ dx. 

i>=n+l 



By division, we find that 

1 + x + x^ ^ ^ ' l^-x + x^ 

IJ.=1 

ioT v ^ 1, where the sum is zero ii v = 1. Therefore, 

, 2 + x A 3^ 1 /3 , \ A^'^A'VS 1 1 
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for n ^ 1, where the double sum vanishes if n = 1. 

To evaluate the last integral in (45), assume that m = 6k < 3"^+^ and let ao,fe, 
0'6k-i,k G ^ denote the coefficients in the polynomial 

6fc-l 

(47) J^^'j^kx' ^{2 + x)Qk{x). 

Noting that ao,fe = (-3)^^=56/? by (44), and using Lemma 1, we have 
(48) 

„1 oo 6fc — 1 oo ^ 

/ {2 + x)Qk{x) J2 x^'^-^dx= J2^j,k Yl 5^ 

•^0 i/=n+l j=0 !/=n+l 

= (-l)'^+y ^Inafl + ^V 

^ 3n-3fc / ^ J ^ 3^ J 

Formulas (45), (46), (48) imply the following base 3 version of Lemma 4. 
Lemma 8. If m < 2t^~^^ is a multiple of 6, say m = 6k, then 

(49) In,m,3 = (-3)^''7 + Ln,ni,3 - ^n,m,3, 

where 

m — l / 

Ln,n.,, = i-3rnlog3+J2 —1^3 1 + 4 

and 



Ar 



V2 3- ^1 3^+1 



1 



3// + 2 



satisfy 



d^Ln,m,3eZlog3 + Zln3[l + ^] +Zln3( 1 + 1- ) +--- + ^ln3(l + ^^^ 



and ci3"An^,^^3 e and t/ie 0-,-,^ integers determined by (47) and (44). 
Lemma 9. We have the bounds 

1 /S^A"" /33\"^ 

(50) -T <^nm3<3- -7 /orm = 3" - 3^0 

and the limit 

(51) lim ^Qg-^'^:3"-3,3 = 3iog3 _ 4log4 = -2.24934057 . . . . 

n— >oo 3" 
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Proof. Following the proof of Lemma 7, with m = 3"^ — 3 we have 

-1 1 /Q3\m, 



(3m + 8)!m!^ (3/4)« /4m\ \^ 1 



and 



(4m + 9)! 4m + 9 V / 40m + 90V44 

/■I o I oo 

" u=n+l 

Jo i + x + x^\ ; 



I/=l 



<[T^) (2 + 7)<3 
This proves the lemma. □ 

Theorem 6. We have the asymptotic formula for Euler's constant 

Proof. Let m = 6/c = 3"^ — 3 for n > 0. Then (—1)^'^ = (— l)**"^ and the formula 
follows from (49) and (51). □ 

Theorem 7. Euler's constant is rational if and only if the fractional part of 
c^3"-f'n,3"-3,3 equals cf3n/„,3n_3^3 for all large n. 

Proof. Multiply (49) by ^3^ and write the result as 

d3nLn,m,3 = d^^In,m,Z + C^3" (^n,m,3 — (" 3)^'^7). 

By the Prime Number Theorem and Lemma 9, for m = 3"^ — 3 and any £ > we have 
cJ3"-fn,m,3 = O ((3"^e^+^/4^)"^) as m — > oo. Thus < ^3^/,^ 3«_3 3 < 1 for all large n. 
Since (isriA^ 3«_3 3 e Z, the theorem follows. □ 

Theorem 8. // the inequality 



{c^3"-f'n,3"-3,3} > 3 ■ ^-^ 



holds infinitely often, then 7 is irrational. In fact, the inequality for a given n > 
implies that no divisor of 3'^ d^n can be a denominator of ^ . 

Proof. The inequalities d^^ ^ d^^N < 8^^ and the upper bound in (50) imply that 

3 3™ 3^~'~''' 
a3"-fn,3"-3,3 < 8 ■ •3-(|4j ^^'(4 

for n ^ 1. Using Lemma 8, this proves the second part of the theorem, which implies 
the first. □ 
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9. Base q asymptotic formulas for Euler's constant 
Denoting by Fq{x), Gq{x) the polynomials 

1 — 

Fg{x) = l+x + x^ + ■■■ + x'^-^ = ^ _ ^ , 

^^^^ o a-F (x) 

Gq{x) = {q-l) + (q-2)x + --- + x^'^ = ^—^^ 

we may write the integral (40) as 

Expanding the binomial and using Ramanujan's formula /o,o,g = 7) we obtain 

+ / Gq{x)J2i-'^y [ jQ'^-'Fqixy-' E ^'""'do;. 
In the first integrand, we find by dividing that 



for f ^ 1, the summation being zero if = 1. Therefore, 

(55) / ?44E-''''-^da. = nfvl-log,VE E " E ^1 

for n ^ 1, where the double summation vanishes if n = 1. 

In the second integrand in (54), let bo^rn,q, • ■ ■ , ^m(q-i)-i,m,Q ^ ^ denote the coeffi- 
cients of the polynomial 

m(q-l)-l m ^ x 

(56) E bi,^,qX^ = Gq{x)Y,i-^yr)Q'^-'Fq{xy-'- 

1=0 j=l ^"'^ 

Note that bo,m,q = {q - l)((g - 1)"" - g"") by (52). If m ^ g"+V(g - 1), then by 
Lemma 1 the last integral in (54) equals 

im(g-l)-l oo m(g-l)-l oo ^ 

(57) / E ^l^rn,qx' E ^'''"'dx= E ^^.-.'J E 

Equations (54)-(57) imply the following generalization of Lemma 4 to base q. 
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Lemma 10. If m ^ q^^^ /{q — 1), then 

where 

m(g-l)-l / 7 \ 

and 

9° If-;/' its v w ^^s/' + a 



/i=l I/ = 2 /[J,: 



satisfy 

dm(q-l)Ln,m,q G Z log g + Z lllg ^1 + -^^ + Z lllg + + 

m{q — 1) — 1 



+ ZlnJ 1 + 



and dq-n.A'^ .j^ q e Z, and the bi^m,q o.re integers determined by (56) and (52). 

To optimize the asymptotics of I'^mq/'f^ ^^e range of Lemma 10, and thus 
obtain the smallest error term in our formula for 7, we choose m ~ (g"' — l)/(g — 1) 
in the remainder of this section. Let fq{x) denote the polynomial 

q-l 

Lemma 11. If q^ 2, there exists Xq between and 1 such that 

(59) fq{xq) = max fq{x). 

Moreover, we have the asymptotic formula 

(60) / fq{x)"'x''-^dx = fq{Xq)'^^''+°^^^^ aS TU ^ OO 

Jo 

and the bounds 



Proof. Since 
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is decreasing for x ^ 0, positive at a; = and negative at a; = 1, it follows that fq{x) 
vanishes at exactly one point the open interval (0, 1). Then (59) holds with 

this choice of Xq, since fq{0) — fq{l) = and fq{x) > for < x < 1. Differentiating 
(62), we substitute fq{xq) = and find that fq{xq) < 0; by Laplace's method (see, 
e.g., [5], Section 4.2) we obtain (60). 
Since 

x + x'^ -\ ha;"-^' 



the inequalities 



q- Fq{x)^{q 



(x • . . . a;'/-i)i/('?-i) ^ x + + ' ^ ^ for ^ X ^ 1 

yield 

(g - 1)(1 - x)a;«(''-^) ^ fq{x) ^ {q - 1)(1 - W^)^^^'^-^). 

Using the relations 

1 r*" 1 

— < max ((1 — v)v^) = r — TT < — for r > 0, 

(r+l)e 0^2/^1^^ ^ (^ + l)r+i re 

we set y = X, r = q{q — 1) and deduce the lower bound in (61), then set y — x^/^, 
r = 2(g — 1) and obtain the upper bound. □ 

The next result with g = 2 is comparable to Lemma 7 with m = 2"' — 1. 

Lemma 12. For q ^ 2, we have the bounds 

q q^ — \ 

(63) < /' „ < -— — with m = — 

V / n,m,q ^2e)"^ q-1 

and 

(64) -1 - iog(, + 1) < '7^!^-;-;;;:-^ < -i - ^- 

Proof. Let m = {q"" - l)/{q - 1), so that a;«"^'-i = a;(^«+i)(«-i). Using (53), (59) 
and Ramanujan's formula -^0,0,(2 — I1 have 

s;/,W"(c, + 7), where c,= [ ^ME>dx, 

Jo ^qyX) 
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and (63) follows using Cg + 7 < — 1) + 7 < g and the upper bound in (61). 
We also have 

(66) /;^,, > j^' ^(g-F,(x))^x(-^+i)(^-i)dx >^^J^'uxrx'^-'dx. 
From (65), (66), (60) we deduce that 



(67) lim ^ W f ) 



and (61) gives (64). □ 

Remark. The exact value of the constant Cq is — ■0(l/g) — logg — 7, as another formula 
of Ramanujan shows (see [2], Proposition 2.1). 

Theorem 9 (Base q Asymptotic Formula for 7). If q 2 and m = {q'^ — l)/{q — I), 
then 

A' — T' 

7 = ■ — ■ — — + (Jm o W^^^'i < dm o < 7 -x 



Proof. Set bm,q = In,mJ(f ^ud usc (58) and (63). □ 

Remark. For = 2, one could use I'm,n,q with m = (g" — l)/(g — 1) to derive 
a rationality criterion and an irrationality test for 7 (but they are already con- 
tained in Theorems 2 and 3 as the case m = 2" — 1, since -/„ nn 1 o — J-ri 2"- 1 

and -^n,2"-i,2 = -^n,2"-i)- For g > 2, this is not possible, the required inequality 
^q^-^n,m,q < ^ being valid, since 

lim i^ = l 

n— >oo g"' 

(from the Prime Number Theorem) and the following observation imply that dqnl'^^^ q 
tends to infinity with n. 

Lemma 13. If q^2>, then 

(68) lim ^""^^'r^M-^yj^-iU > 

n— >oo g" 

Proof. By (64), we have 

^.^ l0S-^n,(q"-l)/(g-l),g ^ 1 ^-^ ^"g -^n,(q"-l)/(g-l),g 
n^oo g" g — 1 n^oo (g*^ — l)/(g — 1) 

1 + log(g + 1) 

> for g ^ 2. 

g- 1 
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Since the lower bound is increasing, and exceeds —1 when q — 4, this proves (68) for 
g > 3. For g = 3, we use (67) to replace the bound by the exact value 

^log/3(a;3) = -0.90997390 ... > -1, 

where = 0.86304075 ... is the positive root of 12 — 7a; — 8a;^, which is the polynomial 
(62) with g = 3. □ 

Remark. Alternatively, for g sufficiently large (68) follows from (67) and the formula 

lim fq{Xq) = ^, 
q^oo Ze 

which in turn follows from the upper bound in (61) and the fact that the limit, as 
g — > oo, of the lower bound in 



max /.(x) ^ Mx)l^ii-i/,y/. = " T^T^StW^) ' " ^) 
is l/(2e). 

10. Concluding remarks 

Any polynomial Pmix) G Z[a;] can be used to generalize integral (18) to base g as 
follows: 



(69) £(_^__l_)p„(.)J^^.,^-.^. 



Integrals (37) and (40) are particular cases of construction (69). The choice of Pm{x), 
the "damping factor," is based on the requirements: 

(i) to ensure the appearance of 7, and no other ^j^g from (39), in the corresponding 
linear forms; 

(ii) to produce sufficiently good asymptotics (upper estimates). 

The choice Pm{x) = (1 — x)'^ gives the best possible answer (38) to feature (ii); that 
is why we use it in Sections 2-4 for g = 2, and Section 8 for g = 3. The following 
lemma and the arguments in Section 7 show that {1 — x)'^ cannot be used for g > 3, 
because of requirement (i). 

Lemma 14. If q > 3, the equality 

(70) [I - x)"^ = C + Q{x) ■ \^ 

1 — X 

is impossible for any m > 0, constant C and polynomial Q{x). 
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Proof. For q > 1, equality (70) implies the system of equalities (1 — e-')"^ = C for 
j = 1,2, . . . , q — 1, where e is a primitive q-th root of unity. In particular, the value 
of |1 — e-' I is the same for all j — 1,2, . . . , q — 1, which is possible only when q — 2 or 
q = 3. D 

One can avoid this difficulty by introducing linear combinations of powers of 1 — a: . 
More precisely, there exist (integer) constants 0^,0^ Cm,i, • ■ ■ , Cm,q-i such that 

c^,o(i - xr + c^,i(i - xr~' + ■■■+ cm,q-i{i - xr-'^+^ = c+ q{x) ■ 

However, we cannot choose the left-hand side as the polynomial Pm{x), because the 
constants grow so much faster than m that requirement (ii) fails. 
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